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Abstract 

Phase transitions in one-dimensional classical fluids are usually ruled out by making appeal to van 
Hove’s theorem. A way to circumvent the conclusions of the theorem is to consider an interparticle 
potential that is everywhere bounded. Such is the case of, e.g., the generalized exponential model 
of index 4 (GEM-4 potential), which in three dimensions gives a reasonable description of the 
effective repulsion between flexible dendrimers in a solution. An extensive Monte Carlo simulation 
of the one-dimensional GEM-4 model [S. Prestipino, Phys. Rev. E 90 , 042306 (2014)] has recently 
provided evidence of an infinite sequence of low-temperature cluster phases, however also suggesting 
that upon pushing the simulation forward what seemed a true transition may eventually prove to 
be only a sharp crossover. We hereby investigate this problem theoretically, by three different and 
increasingly sophisticated approaches (i.e., a mean-field theory, the transfer matrix of a lattice 
model of clusters, and the exact treatment of a system of point clusters in the continuum), to 
conclude that the alleged transitions of the one-dimensional GEM4 system are likely just crossovers. 
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I. INTRODUCTION 


Compared to three dimensions, the statistical mechanics of one-dimensional (ID) fluids is 
both simpler and harder. Simpler, because a general result due to van Hove [1] excludes the 
possibility of distinct phases in one-dimensional, homogeneous, and one-component classical 
fluids of hard particles with short-range interactions. Also more complicated, however, 
because whenever any of the hypotheses behind van Hove’s theorem is not met, it may be 
extremely difficult to say whether or not the given ID interaction would admit a genuine 
phase transition (see, e.g., Refs. M). 

In this respect, the class of potentials known as generalized-exponential models (GEM) 
is emblematic. These are classical fluids of particles repelling each other through a potential 
of the form u{r) = eexp{ —(r/a)"} (GEM-a potential), where e,a > 0 are arbitrary energy 
and length units, r is the interparticle distance, and a is a real positive index. The GEM 
potentials stay hnite at the origin, hence particles may fully overlap, in strike contrast to real 
atoms and molecules. Glearly, bounded potentials only make sense as models of the effective 
pair interaction between large molecular aggregates that can intertwine even to the extent 
that their centers of mass exactly coincide. In three dimensions bounded interactions give 
rise to structural patterns which are unknown to rare-gas fluids, notably reentrant melting 
and waterlike anomalies, which have been the subject of numerous investigations in the last 
few decades (see, e.g., Refs. to name but a few). A particularly curious feature of 

bounded two-body repulsions is the possibility to generate stable cluster crystals [15], i.e., 
crystalline phases involving a mean number of particles per lattice site larger than 1, which 
have been extensively studied in the recent past [T6lj2l] . 

When conhned in one dimension, penetrable particles may pass one through the other, 
thus making it possible to change the particle ordering within the line. While there exists 
a general method of computing the exact thermodynamics of a specihc class of ID systems 
of impenetrable particles [25H2Z] , nothing can a priori be said about the phase behavior of 
freely overlapping particles. Hence, either a brute-force numerical-simulation approach is 
taken or some other plausibility arguments must be used to probe the existence of separate 
phases in such ID systems. 

In Ref. [8], the ID GEM-2 (that is, the one-dimensional counterpart of the Gaussian- 
core model [281 - I3T] ) was studied by Monte Garlo (MG) simulation. No phase transition 
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was found, even at zero temperature (T = 0), notwithstanding the system structure at 
low temperature is crystal-like and clear signs of dimerization were found at high pressure. 
For the ID GEM-4 the situation is markedly different [ 32 ] • True phase transitions exist at 
T = 0 between states allowing for an increasing number of particles in the same lattice 
site. In simulation, these states gave rise to clear-cut cluster phases for T > 0, with sharp 
boundaries which manifest themselves through hysteresis loops in the number density as a 
function of pressure. In fact, the coexistence loci are so sharp that they could be traced 
by free-energy calculations. However, when MC simulations were made much longer (up 
to tens of million sweeps), hysteresis loops became narrower and the alleged critical points 
shifted to low temperatures, thus raising doubts on the true nature of the crossover from 
one phase to the other. Therefore, it is still unclear whether proper phase transitions truly 
exist or not in the ID GEM-4. 

In the present paper, we address this problem from the different perspective of theoretical 
modeling. Upon building up successive, increasingly accurate approximations to the ID 
GEM-a system for a > 4, we compute the density as a function of pressure for hxed 
temperature, seeking for jumps or other singularities. At the end, we get a more dehnite 
idea about the nature of phase crossovers in GEM fluids. Even though not having the same 
degree of certainty as a mathematical proof, nonetheless the present work sheds fresh light 
on a still open question. 

This paper is organized as follows. After introducing a coarse-grained formulation of the 
ID GEM-a Hamiltonian in terms of interacting clusters on a lattice, we compute the phase 
boundaries of the model hrst in the mean-held approximation (Section H), and then by the 
exact transfer-matrix method (Section HI). Going back to the continuum, in Section IV we 
analyze the exact behavior of a system of point clusters moving on a line and discuss its 
relevance for the phase behavior of the original GEM system. Some concluding remarks are 
given in the last section. 


II. MEAN-FIELD THEORY 

The statistical mechanics of the ID GEM-a system is not amenable to an exact an¬ 
alytical treatment. However, some progress can be made by analyzing a related system 
which allows to simplify the counting of states. In this regard, a crucial observation from 
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the simulation |32] is that the low-temperature structure of the GEM fluid away from the 
alleged “transition points” is not dissimilar from an ordered cluster crystal, with average 
site occupancy and lattice spacing depending on temperature T and number density p in a 
smooth way. At a phase transition, the average occupancy increases by one and the density 
undergoes a sharp variation (however, close to a transition point clusters do not necessar¬ 
ily contain the same number of particles and the separation between two adjacent clusters 
may then vary along the line). Given the above evidence, we are led to represent the low- 
temperature behavior of the ID GEM-a system of N particles in a volume L by means of 
the coarse-grained lattice model of Hamiltonian: 

^ Afc Vc 

H = -e'^ni{ni - 1) + rij u {m tnn) rij+rn , (2.1) 

2=1 2=1 m >0 

where Nc is the (for the moment arbitrary) number of point-like clusters (one for each lattice 
site), tnn = L/Nc is the lattice spacing, he., the distance between nearest-neighbor (NN) 
clusters, and = 1, 2, 3,... is the size of the cluster residing at i, he., the number of particles 
located at the hxed lattice positions 


Xi = {i- l)rNN (f = 1,..., A^c) • 


( 2 . 2 ) 


On the right-hand side of Eq. (2.1), the hrst term is the intra-cluster energy while the 
second term represents the total interaction energy of particles belonging to different clusters 
(periodic boundary conditions are implied). Observe that the inter-cluster potential in 
Eq. ( |2.1 ) is not restricted to nearest-neighbors, but involves arbitrarily distant clusters. A 
model similar to (2.1) has already been considered by Wilding and Sollich |33] as a tool to 


investigate the T = 0 transitions of the three-dimensional GEM-a system. In this section 
we extend their mean-held treatment to T > 0 for the ID case, while deferring the exact 


analysis of model (2.1) to the next section. 

In order to extract the thermodynamic properties of ( |2.1 ) we hrst use the variational 
principle of statistical mechanics. Since handling the rii = N constraint is hard, we 
shall work in the grand-canonical ensemble where T, L, and the chemical potential /i are 
hxed, while each Uj can assume every value between 1 and oo. In this ensemble the Gibbs- 
Bogoliubov inequality prescribes that the grand potential D is bounded from above by a 
functional D*, 


D < D*[7rl = Tr 


7r(if -f — liiTT — ^N) 


= E* - TS* - nN *, 


(2.3) 
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where vr is any trial probability density in state space and {3 = hi* reaches its 

minimnm for tt = Pg.c. (he., the eqnilibrinm grand-canonical density), where it takes the 
valne hi. Once a parametric form of tt is chosen, minimizing O* [vr] with respect to parameters 
is the way to make the best possible estimate of O. 

A system with no specihc constraint on the total nnmber of lattice sites entails a com¬ 
plication. In fact, after compnting O (he., the least npper bonnd to O) for each Nc, we still 
need to optimize with respect to N^. This is done according to the criterion stated by Swope 
and Andersen |21], that is by minimizing O with respect to = L/N^, or, eqnivalently, 
with respect to the Ac-dependent density d = N*/L [53] . 

The simplest choice is a 7i{ni,n2, ■ ■ ■, of nncorrelated n* variables. Following Ref. [33] . 
we choose n (he., the mean site occnpancy) as the only parameter in vr. At low temperatnre, 
we are allowed to assnme that only two n* valnes at a time may occnr in the system, namely 
= int(n) and -1- 1. Hence we have: 

Vc 

7r({n};rr) = Jj7ri(ni;n) , with 7ri(n; n) = (1 - -h . (2.4) 

i=l 

The valne of A is determined by the condition = Tr(7mj) = n. Since 

{ni)-K = ^ni7ri(?7.i;n) x ^7ri(n2;n) x ■ ■ ■ = (1 - A)?7.< -h An> , (2.5) 

ni 712 

we get A = n — n^, which also yields: 

{^ 1)77 = (1 ~ ^)’^< + =n‘^ + A(1 — A). (2.6) 

From the very dehnition of N* it follows that N* = NcU, hence 

N^ = Ld/n. (2.7) 


Pntting all things together we obtain: 


Ld 


A* = — < -I- A(1 — A) — nl e -I- 2n^ 

2 n 

Ld 


n\ (2n 


-rs* = ksT— [(1 - A)lii(l - A) + Ain A] ; 
—fiN* = —fiLd , 


( 2 . 8 ) 


and the system pressnre P is then given by: 

P(T, /i) = — — minr2((i; T, p) = —min min Q*{n, d] T, 
L d Ld 


(2.9) 
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Any jump in the equilibrium number density p = argmin^ T,/i)} as a function of p 

at constant T will be the hallmark of a hrst-order transition. 

In Figs. 1 and 2 we have plotted our mean-held results for the ID GEM-4 system relative 
to six isothermal paths, T = 0, 0.01,..., 0.05 (from now on, all quantities will be given in 
the units set by /crA, and a). For these temperatures we hrst computed the best value 
of n and the corresponding value F oi F* = E* — TS* as a function of d (top panel of 
Fig. 1). The optimal n undergoes a steady increase with d, interrupted at regular intervals 
by plateaus at integer heights. Within the same d ranges where n stays constant, F is 
convex upward. The bottom panel of Fig. 1 shows the lattice spacing r^N = L/Nc = n/d 
as a function of d; tnn exhibits a zig-zag pattern which gradually dampens as d grows. 
Observe that the value of tnn never falls below 1.1, meaning that the distance between two 
next-nearest-neighbor clusters is always larger than 2.2. Since exp{—2.2^^} < 10“^°, this 
means that the GEM-4 interactions are effectively limited to nearest neighbors only (the 
same applies for any a > 4). In Fig. 2 we show p (top panel) and p (bottom panel) as a 
function of P for the same temperatures considered in Fig. 1. The density jumps at the 
lowest temperatures indicate a sequence of hrst-order transitions between cluster phases of 
increasing occupancy. Above T 0.03 all discontinuities vanish and the transitions are 
replaced by smooth crossovers. 

Summing up, our mean-held theory asserts the existence of low-temperature phase tran¬ 
sitions in the ID GEM-4 system. Gompared to the work by Wilding and Solhch[33], two 
important improvements are: 1) our theory is directly derived from the variational principle 
of statistical mechanics and, as such, it can be systematically improved; 2) at variance with 
the theory in Ref. [33j, to which ours reduces for T = 0, the present mean-held theory also 
applies for T > 0. 


III. TRANSFER-MATRIX TREATMENT 


Since the results of an approximate theory may not be conclusive, the accuracy of our 
mean-held theory for ID GEM-a systems will now be tested against an exact treatment of 


model (2.1). We start by rewriting the latter Hamiltonian as: 


Nc 


Vc 


H = - 1) + ^ t((rNN)n,ni+i, 


(3.1) 


2 = 1 


2 = 1 
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which assumes zero interaction beyond nearest neighbors (we shall check a posteriori that, 
for a > 4, the omitted terms are actually ineffective at all temperatures of interest). Again, 


the clusters are fixed at the lattice positions given by Eq. (2.2), and tnn = L/N^ is an 
adjustable parameter to be eventually optimized by requiring the system pressure to be 
maximum. Each Ac-uple of rij values (for i = l,...,A"c) represents an individual state 


of the model (3.1), since the particles building up the clusters are indistinguishable. The 


canonical partition function would be the sum of the Boltzmann factors over all A^c-uples 
(rii,... ,nArJ of positive integers obeying the constraint Therefore, as in the 

previous section, the calculation is most easily performed in the grand-canonical ensemble, 
since the partition function of the latter involves an unconstrained sum over ni,..., 

Using the transfer-matrix method, the grand potential (under periodic boundary condi¬ 
tions) reads 

-/3fi(T,L,/r;rNN)=lnTr{e^^^e-^^}=ln 

= ln5^(T'^^);;; ~iV,lnA^ax, (3.2) 


ni 


where the last step holds for large A), (and L) values. In the above equation, A^ax is the 
largest eigenvalue of the (symmetric) transfer matrix 

+ n'"^ — n — n' 


T", = exp -(5e 


-|- /5/i(n -|- n')/2 — (3u{r-fi-i<i)nn 


(3.3) 


It is worth stressing that the thermodynamic limit implies taking not only L —)• oo but also 
Nc ^ oo, m such a way that L/Nc = tnn remains finite. Also observe that T is an infinite 
matrix, since both n and n' can vary between 1 and oo. In practice, if the system density 


is not too high, we can assume n and n' in Eq. (3.3) to vary between, say, 1 and 10 (T thus 
becomes a 10 x 10 matrix). 

The calculation of A^ax was carried out numerically by the power method [HB]. First 
the system pressure for every Tnn was computed through r2(T, L, /r; Tnn) = —'P{T, /r; Tnn)-^, 
yielding in the thermodynamic limit; 

hr Amax 


V{T, jj,] r nn) — 


l3r NN 


(3.4) 


Such a quantity was then maximized to get the equilibrium pressure: 


F(T, fi) = max F(T, /r; Tnn) • 

f’NN 


(3.5) 
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A first-order phase transition occurs whenever the P(T, p;rNN) vs. tnn prohle exhibits two 
different maxima, whose relative stability changes when crossing a specihc /i value. This 
implies a jump discontinuity in the optimal value of r^N, which in turn produces a jump in 
the density as a function of /i or P. 

We hrst computed P for T = 0.01,..., 0.05, hnding values very close to those given by 
mean-held theory (pressure differences in reduced units were typically smaller than 10“^); 
only for the smallest T the numerical calculation of the pressure could not be carried out 
above a certain p, because of the huge values of some T components. We verihed that, 
even including second-neighbor terms in the Hamiltonian, the pressure does not change 
appreciably for a > 4. The reason is the same mentioned previously: the optimal tnn 
values are never smaller than l.la, hence the strength of the second-neighbor repulsions 
is nearly zero. In Fig. 3 we show our exact transfer-matrix results for the lattice model 


(3.1). Besides tnn, the top panel of Fig. 3 also reports the connected correlation between 
NN clusters, {nin 2 )c = — {ni){n 2 ), and the correlation length ^ dehned by 

? = I |T ,, I . (3.6) 

In |Amax/A2I 

where A 2 is the second (next to the leading) eigenvalue of T in descending order. While 
A 2 was computed by Wielandt deflation 1371, the NN correlation was determined by the 
formula: 


(^ 177 , 2 ) = 


A. 


10 

E 


pniTn,n2 




(3.7) 


ni,n2=l 

where the M matrix has the T eigenvectors for columns (in descending order of eigenvalue). 
The absolute value of {nin 2 )c, which is very small for all pressures, may explain the good 
quality of our mean-held theory. The bottom panel of Fig. 3 reports the average occupancy 
(ni), the average squared occupancy (n^), and the density p, respectively computed by the 
formulas: 

10 10 

(w) = ’ (^ 1 ) = X] and p = (ni)/rNN • (3.8) 


ni=l 


ni=l 


Overall, the density behavior closely resembles that found in mean-held theory. Indeed, we 
see again discontinuous jumps at the lowest temperatures (T = 0, 0.01,0.02), indicating true 
hrst-order phase transitions. At higher temperatures, the discontinuous jumps transform 
into sharp, but continuous, density changes in a very narrow range of pressures (similar 
trends were named pseudo-transitions in a diherent ID model by P§kalski et al. 0 ). 





The correlation length shows a distinct maximum near each low-temperature transition 
point. As T grows from 0, these maxima hrst increase and then, past T = 0.03, start to 
decrease, suggesting that each coexistence line ends with a critical point and that ^ diverges 
at a critical temperature ~ 0.03. Since the transfer matrix T is strictly positive for any T, 
one might argue that its maximum eigenvalue is necessarily simple (by Perron-Frobenius 
theorem |38]), and therefore it would be safely excluded that ^ may ever diverge. However, 
as already noted, the true transfer matrix is infinite and therefore the use of the Perron- 
Frobenius theorem is not legitimate [35]. Hence, the question of the behavior of ^ near 
T = 0.03 can only be answered numerically (see below). 

Upon looking more deeply into the density and correlation-length profiles close to every 
zero-temperature transition point, we observed a very complex behavior. In Figs. 4 and 5 
we focus on the transition from the 1-cluster phase (1) to the 2-cluster phase (2). The 
density jump at the transition has an irregular T behavior. Upon heating from T = 0, 
the 1-2 coexistence line “bifurcates” at T < 0.028. Indeed, two distinct density jumps are 
seen for T = 0.028 (red curves in Figs. 4 and 5), and this entails the existence of a third 
phase, intermediate between 1 and 2, characterized by values of n < 1.5. However, this 
extra-phase has a very precarious existence, since it has already disappeared for T = 0.029. 
Eventually, above T = 0.033 we find no density jump but a continuous density curve, and 
the 1-2 coexistence line terminates. Figure 5 shows that the correlation length never diverges 
and that the Widom line {i.e., the locus of regular ^ maxima), rather than departing from 
the ending point of the coexistence line, actually blossoms from the coexistence line near 
T = 0.028. We hnally notice that the behavior of p and ^ near the 2-3 and 3-4 transitions 
is perfectly analogous to that described for the 1-2 transition. 

Figures 6 and 7 report the overall phase diagram on the P-T and p-T plane, respectively. 
On each coexistence line there is a bifurcation point and the Widom line apparently emerges 
from it. The mean-field coexistence lines are also shown. They run extremely close to the 
exact loci up to the bifurcation points, while progressively departing from them for higher 
temperatures. 

Finally, we looked at what happens when increasing a from 4. We found that the discon¬ 
tinuous density jumps are only present in a smaller and smaller range of low temperatures 
until they altogether disappear for all T > 0 when a —)■ oo. In this limit, (ni) coincides with 
p (see Fig. 8). Hence, the ID penetrable-sphere model (PSM, corresponding to the a —)■ cxd 
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limit of GEM-a) would show no phase transition at T > 0. 


IV. A CONTINUUM MODEL OF INTERACTING CLUSTERS 


We hereafter show that in Hamiltonian (3.1) the lattice constraint (2.2) can be weakened, 


making the distance between clusters partly arbitrary, and the model still remains exactly 
solvable. 

In the previous section, we considered a model whose states are represented by W-uples 
of positive n* integers, dehned at hxed lattice positions Xj. Slightly abusing the notation. 


a model state was any vector (ni,a;i,..., xatJ, with Xi dehned by Eq. (2.2). We now 
release the constraint on the cluster positions and let Xi be any real number in the interval 
[0, L]. We focus on the specihc order xi < X 2 < ■ ■ ■ < xn^, which is one out of the NJ possible 
orderings of coordinates. Therefore, a microstate of the by-now continuum model is going 
to be represented by a 2V(.-uple (ni, xi... ,nN^,XNj, with 0 < Xi < X 2 < ■ ■ ■ < x^^ < L. 

Since the clusters are now free to move, the new Hamiltonian is written as: 


Nc 


Nc-1 


p2 ^ 

hf = ^ + o'" X] “ 1) + Uaixi+i - Xi)nini+i , 

ZiTTli L, . 

2=1 2=1 2=1 


(4.1) 


where p* and m* = nim are the momentum and mass of the i-th cluster, respectively (m 
being the mass of an individual particle). In the interaction term, which is again limited to 
nearest-neighbors, the potential Ua differs from u for the addition of a suitable hard core of 
diameter a. We need such a modihcation to ensure that, for the given u potential, a cluster 
effectively interacts only with its hrst neighbors; under this assumption, the calculation 
of the partition function becomes feasible by elementary methods, as we shall see below. 


Admitting only NN interactions in Eq. (4.1) actually amounts to cut the u interaction at 


r = 2a. For the GEM-a potential with a > 4, this is clearly allowed for a > 1.1, but 
we can take an even smaller a if, in the system conhgurations with highest Boltzmann’s 
weight under the given thermodynamic conditions, the number of second-neighbor clusters 
at distance 2a is statistically irrelevant (as occurs in GEM-a systems, with a = 4 or higher, 
for low temperatures). In other words, replacing u with Ua would not alter the statistical 


properties of the GEM-a model, in so far as the most relevant conhgurations of model (4.1) 
are those where the distance between NN clusters is around 1.2 1.5. 
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For reasons which will be apparent below, it is far easier to compnte the partition fnnc- 


tion of model (4.1) in an ensemble where the independent variables are T, P,/r, Nc- Observe 


that T, P, /i do not exhaust the set of intensive variables for our system since fic, the chem¬ 
ical potential conjugate to Nc, is also intensive. Denoting 3^ the partition function in the 
T, P, /i, Nc representation, the appropriate thermodynamic potential is 


6(T, P, fi, iV,) = -keT In = fie(r, P, p.)N ,. 


(4.2) 


However, thermodynamic equilibrium requires that fic{T,P,fi) must be identically zero [201 
[Ml 140114T] and Q then vanishes exactly. It is worth stressing that the equilibrium condition 


Pc{T,P,lJ.) = 0 


(4.3) 


represents itself an implicit functional relation between T, P and p, and we will exploit it to 
express, e.g., p as a function of T and P. 

Statistical mechanics allows to compute 3^ from the Hamiltonian. All details of such a 
calculation are given in the Appendix. The main result is that, similarly to the treatment 
of Section HI, also y{T, P, /i, Nc) can be expressed in terms of a transfer matrix T (defined 


at Eq(A.4)). In the thermodynamic limit, y reduces to 

3 ^ = . 


(4.4) 


LoA{(3P)^ ’ 

where Amax is the maximum eigenvalue of T (and the other, less important, symbols are 


dehned in the Appendix). Substitution of this result into Eq. (4.2) then yields 


Q{T,P,ll,Nc) = -/CBPA^clnAmax + 0(1 ), 


(4.6) 


and from Eq. (4.3) we can conclude that 


Xm^JT,P,f^) ^ 1 . 


(4.6) 


We have thus demonstrated that the equilibrium condition based upon the chemical poten¬ 
tial, He = 0, is perfectly equivalent to the condition Amax = 1 for the largest eigenvalue of 
the transfer matrix T. 


Before going on, it is worth remarking that the same Eq. (4.6) can also be obtained by 


a different route, which is reminiscent of the treatment of ID multi-component mixtures 
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by Longuet-Higgins [13H15] . In this case, let us start from the grand-canonical partition 
function S, dehned as 


j(T, h, p,/ic) ^ ^ 




E 


./3m Ei 


/Vc = l 


ni,...,njVc 




d^d^^xe"^^. (4.7) 


Assuming Pc = 0 from the outset, the partition function actually reads: 

1 






Nc=l ni,...,njvc 

For any number p, let us consider the integral: 


Jln. 


d%d^^xe-^^. 


(4.8) 


I = — dLe-f^P^E{T,L,p,0) 
^0 Jo 


(4.9) 


Denoting P{T, p) the yet unknown system pressure in terms of T and p, and observing that 
S = exp{/3P(T,/i)L}, it follows for any p > P{T, p) that: 


I = — dL e-/3(p--P(EM))i- = 


^0 Jo 


^Lo[p-PiT,p)] 


(4.10) 


On the other hand, by the same steps leading to Eq. (4.4) the integral in Eq. (4.10) can also 
be more cumbersomely written as: 

.,2 


^=E 




2 \ Vc-l 


E 


wz 


(4.11) 


Nc- 


LoA(/?p)2 ^ LoAifJpy tt 1 - An ’ 

{An} being the set of (real) T eigenvalues listed in decreasing order (notice that all An are in 


modulus less than 1 for p > P(T, p), since otherwise one of the series in Eq. (4.11) would not 


converge, contrary to what implied by Eq. (4.10)). Comparing the last two equations, we 
see that when p attains the value P{T,p) the largest eigenvalue Amax = Ai(T, p, p) becomes 


equal to 1, meaning that P = P{T,p) is implicitly dehned by Eq. (4.6). 

Clearly, the just stated condition would make the thermal properties of the model depend 
on the exact value oi X = h/\/27rmea'^, a quantity which sets the value of A^Jk■QT/e relative 


to a (see Eq. (A.4)). We suspect that this annoying feature of model (4.1) is the price to pay 
for clusters to be considered as point-like particles. In the following, we take the pragmatic 
view to choose the X value which produces the best possible agreement with existing MC 
simulation data for the ID GEM-4 system. 

For a = 4, we computed p and p = {dP/dp)T as a function of P in the interval 0.75- 
0.9, corresponding to the 1-2 transition region, for a number of temperatures between 0.002 
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and 0.025. As a preliminary check, we verified that the p valnes are relatively insensitive 
to a: with either a = 0.5 or a = 1 the resnlts were practically the same, snggesting that 
the distance between NN clusters never falls below 1 in the configurations with the highest 
statistical weight. This result is very important, since it is then immaterial which potential, 
either the original u or Ua (for a > 0.5 and NN interactions only), is used for describing the 
properties of the GEM-4 [46j . 


Our results are plotted in Fig. 9 (top panel), for X = 0.01 and a = 1. Upon comparing 
the present results with MC density data for T = 0.02 [32], we see that the improvement of 


the continuum model (4.1) over the lattice model of Sect. Ill is dramatic (a few p curves for 
the lattice model were reproduced for convenience in the lower panel of Fig. 9). In fact, there 
is a slight dependence of the location of the 1-2 transition on X (all density curves are rigidly 
shifted to the right as X is reduced), but this effect is smaller for the lower temperatures 
and actually negligible near T = 0. However, the crucial point is that, compared to the 
lattice model of Sect. Ill, the scenario set by the continuum model is qualitatively different: 
the 1-2 transition, which for low T is discontinuous in the lattice model, is now turned into 
a smooth crossover. This conclusion is consistent with heuristic arguments forbidding any 
sort of spontaneous symmetry breaking in one dimension jlT] HHj. One might observe that 
this outcome does not come as a surprise, and could have been anticipated from the very 
beginning, since endowing the interaction potential with a hard core restores the validity 
of van Hove’s theorem. Although correct, however, in no way this reasoning diminishes 
the virtue of a model whose surprising effectiveness, along with the insensitivity of density 
curves to a, make us believe that it is indeed capturing the very essence of the ID GEM-4 
system. 


The noteworthy efficacy of the continuum model would cast a shadow on the relevance 
of the lattice model for the behavior of ID GEM-4. In other words, the phase transitions 
discussed in Sect. H and HI now appear to be an artifice of assuming a strict crystal struc¬ 
ture for an originally continuous system of low dimensionality. Nor the simulation results 
of Ref. [32] can be taken in support of first-order transition behavior. In fact, the den¬ 
sity hysteresis loops observed in the simulations appeared to progressively narrow as the 
MG trajectories became longer and longer, thus rendering those numerical results actually 
inconclusive. 
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V. CONCLUSIONS 


Particles interacting through pairwise repulsions bounded at the origin are known to 
exhibit unconventional thermodynamic and structural behavior, such as reentrant melting, 
waterlike anomalies, and cluster crystals. Which features are observed in a given case 
crucially depends on the range of the interaction and the space dimensionality. 

In three dimensions, the generalized exponential model of index a (GEM-a) shows reen¬ 
trant melting for a = 2, whereas it gives rise to a sequence of stable cluster-crystal phases 
for any a > 2. In this paper, we have studied the phase behavior of the one-dimensional 
(ID) GEM-a systems by theoretical means, focusing on the a = 4 case where numerical- 
simulation results are available |32]- MG simulations gave contrasting indications: on the one 
hand, the low-temperature hysteretic behavior of the number density points to the existence 
of low-temperature first-order transitions between cluster phases; on the other hand, the 
systematic narrowing of hysteresis loops, as simulations were made longer and longer, sug¬ 
gests that the apparently singular behavior of the ID GEM-4 system is actually a finite-size, 
hnite-time artifact. 

To clarify the situation, we introduced both a lattice and a continuum cluster model 
aimed at specifically representing the low-temperature behavior of the original GEM system. 
While the lattice system exhibited a sequence of low-temperature first-order transitions 
between cluster phases of increasing occupancy, as evidenced by the discontinuous jumps 
of the number density as a function of pressure, the more accurate model defined on the 
continuum only showed a series of smooth crossovers between cluster “phases” devoid of 
separate individuality. When comparing the results for the continuum cluster model with 
MG simulation data for a = 4 we found a rather good agreement, thus indicating that strict 
phase boundaries are likely absent in the ID GEM-4 system, and even more so in the ID 
GEM-a with a > 4. 
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Appendix A: Calculation of y{T, P, fi, Nc) 


The kinetic part of y is immediately calculated after noting that 


1 / dp , exp |- MU ^ 


(A.l) 


where A = h/\/2wmk^ is the thermal wavelength of an individual particle. The partition 
function then becomes: 


y{T,p,fi,N,) = 


y Lo Jo A., 






E 

ni,...,njvc 

/-oo 


g/3MEi ’T’ig-AA2) E* ni{ni-l)^{1/2) J2i Inni 


""LoJo 


dXi ■ ■ ■ dXN, g-/3Ei«a(x,+i-xPn,n,+i ^^^ 2 ) 


where Lq is an arbitrary length. 


In order to compute the integrals in Eq. (A.2), we shall proceed d la Takahashi [25], writ¬ 


ing the inner integral as a chain of convolutions whose Laplace transform is easily computed 
(that is why working at hxed P is more convenient). The result is (cf. appendix A of 
Ref. Hg): 


dXi • • • dXN, ^-0j:,u4x,+^-xPn,ni+, 


^ rdLe-^^^^ 

Lo Jo a.. 

AqA ^ Jo J xi<...<x{f^ 

1 A^c-l fOO 1 

_ JJ' / _ p-/3Ps^-/9tta(s)nini+i 


LoA(/?P) 2 7o A 


(A.3) 


Notice that Eq. (A.3) is only valid under the assumption of open boundary conditions. Upon 


defining the transfer matrix T as the following symmetric matrix, 

rjpn _ /3/i(n+n')/2 -/3e(n2+n'2-n-n')/4 (l/4)(lnn+lnn') f -fSPs -I3u{s)nn' 

11„,-e e e / T 


(A.4) 


Eq. (A.2) is eventually rewritten as: 
1 


3^ = 


AoA(/?P)2 


g/3/4ni/2g-^eni(ni-l)/4g(l/4) lnnirjpnirjpn2 , , , qp^jVe-i 


ni,...,njVc 


^ g/3MniVc/2g-^enjVe(niVe-l)/4g(l/4)lnn^^ _ 


LoAifdpy 


vTt^--U 


(A.5) 
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where v is the iVc-dimensional column vector with components 


y. — g^M>^i/2g-^eni(ni-l)/4g(l/4)lnnj 


(A.6) 


Equation (A.5) is an exact result, which involves Nc in a complicate manner. However, in 


the Ac —)■ oo limit Eq. (A.5) is radically simplihed, as we are going to show soon (observe 
that, when Ac —)■ cxd also A = '^)- 


Let M be the orthogonal matrix having the T eigenvectors for columns. Then 
vTjW.-iv ^ with w = R“‘v and r = R“'TR. 


(A.7) 


Called Amax the maximum eigenvalue of T, in the thermodynamic limit the expression of y 
simplihes to 


3^ = 


w 


1 \iVc-l 


LoA(/3P) 


2 ■^max > 


(A.8) 


which is the same as Eq. (4.4). 
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6 



FIG. 1. (Color online). ID GEM-4, mean-field results for six temperatures, T = 0,0.01,..., 0.05, 
marked as black, blue, cyan, green, magenta, and red lines, respectively. We here show results 
deriving from the optimization of n only (top: best n and F/L values; bottom: lattice spacing). 
Each inset shows a magnification of the d interval from 1.3 to 1.7 (top: n; bottom: lattice spacing). 
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FIG. 2. (Color online). ID GEM-4, mean-field results for the same temperatures considered in 
Fig. 1. We here plot the density p (top) and the chemical potential p (bottom) as a function of 
the pressure P, as resulting from the optimization of both n and d (only the first three transitions 
are shown). In each inset, we show a magnification of the P interval from 2.3 to 2.7, which makes 
us better appreciate the effect of heating the system: while the transition from the 2-cluster to 
the 3-cluster phase is still first-order for T = 0.03, it has by now become a sharp (but smooth) 
crossover for T = 0.04. 
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FIG. 3. ID GEM-4, transfer-matrix results for various temperatures (T = 0,0.01,... ,0.05, same 
colors as in Figs. 1 and 2). In the top panel, the average NN distance, the correlation length, 
and the connected NN correlation are plotted as a function of pressure. In the bottom panel, the 
average cluster size (ni), the average squared occupancy, and the density are shown. 
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FIG. 4. ID GEM-4, transfer-matrix results for various temperatures (T = 

0.01, 0.02, 0.025,0.026,..., 0.034, 0.035,0.04,0.05; for the sake of clarity, the line colors are 
alternatively black and blue, except for T = 0.028 ~ red curves - where the density shows two 
distinct jumps at two nearby pressures): average cluster size (top panel) and number density 
(bottom panel) near the transition from the 1-cluster to the 2-cluster phase. In the inset, we show 
a magnification of the density for T = 0.028, which clearly signals the existence of two distinct 

steps (the raw data are the dots, joined by straight-line segments to help the eye). 
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FIG. 5. ID GEM-4, transfer-matrix results for various temperatures (the same as in Fig. 4): 
correlation length ^ near the transition from the 1-cluster to the 2-cluster phase. 
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FIG. 6. ID GEM-4 phase diagram on the P-T plane: Exact transfer-matrix results (open blue 
circles and squares) vs. mean-field results (black dots). The exact Widom line is also shown (dotted 
red line). 
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FIG. 7. ID GEM-4 phase diagram on the p-T plane: Exact transfer-matrix results (open blue 
circles and squares) vs. mean-field results (black dots). 
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FIG. 8. ID GEM-oo (PSM), transfer-matrix results for various temperatures (T = 0,0.01,..., 0.05, 
same colors as in Figs. 1, 2, and 3): average cluster size (ni) (which coincides with the density) 
and average squared size (nf). In the a —)• oo limit, the discontinuous jumps of the density at the 
lowest temperatures have disappeared, leaving a sharp but continuous crossover at all non-zero 
temperatures. 
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FIG. 9. ID GEM-4, comparison between /r and p results for the continuous model (Section IV) and 
for the lattice model (Section III) in the 1-2 transition region (for X = h/VPirmaP = 0.01 and 
a = 1). Top: continuous model (T = 0.002,0.003,0.004,0.005, black; T = 0.01, blue; T = 0.015, 
cyan; T = 0.02, green; T = 0.025, magenta; T = 0.03, red). Observe that even for such a low 
temperature as T = 0.002 the density varies smoothly across the “transition”. In the same panel, 
the black dots are T = 0.02 density data from Ref. [32]. Bottom: lattice model (T = 0.01, blue; 
T = 0.02, green; T = 0.025, magenta; T = 0.03, red). 
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